Modeling interaction of trading volume in financial dynamics 
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Abstract 

A dynamic herding model with interactions of trading volumes is introduced. At time t, an agent trades with a 
probability, which depends on the ratio of the total trading volume at time t - 1 to its own trading volume at its 
last trade. The price return is determined by the volume imbalance and number of trades. The model successfully 
reproduces the power-law distributions of the trading volume, number of trades and price return, and their relations. 
Moreover, the generated time series are long-range correlated. We demonstrate that the results are rather robust, and 
do not depend on the particular form of the trading probability. 
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1. Introduction 

In statistical analysis of financial markets, much attention has been drawn to the study of the stock prices |Q1 0, 
S H H S S H Isl E EH EH EHl- Denoting p(t) as the price of a given stock or financial index, the price return 
r(f) is defined as the change of the logarithmic price in a time interval At, i.e., r(t) = \np{t) - In p{t - AfV A 
power-law tail with an exponent g r = 3.0 is found in the cumulative probability distribution of price returns SHIH], 
which indicates that the large price fluctuation is more common than one might naively expect. Another important 
statistical property basically observed in financial markets is the long-range correlation of the volatility which is 
simply defined as the magnitude of the price return li. Many efforts have been devoted to the understanding of 
the financial markets along this direction, with both phenomenological analysis and microscopic multi-agent models 

d! S B S H EI M EH m M M EH lil . 



Recent empirical studies show that the trading volume is highly correlated with the price return and volatility 
IH |23l |3 H, and this confirms the famous saying that it takes trading volume to move stock pric es. A positive 



linear correlation is revealed based on the data analysis at time scales large than one minute 1 26, 27|,|28|]. For the high 



frequency data at microscopic transaction level, the volume-return relation follows a scaling behavior. Lillo et al 



found a mater curve with scaling form using the Trade and Quote database of US stocks [29], and the scaling function 
is found to be a power-law form for large volumes. Lim and Zhou found similar scaling behavior in Australian and 
Chinese stocks 1 3(1 31]. Due to the significant importance of the trading volume, its statistical properties is worthy of 



being carefully analyzed. 

Interestingly the cumulative probability distributions of the trading volume as well as the number of trades also 
show power-law behaviors, and the corresponding exponents are reported to be £y = 15 and = 3.4 respectively 
(0, HI |25t]. To understand the power-law distributions of the trading volume and number of trades, and their possi- 
ble relations with the distribution of price returns, an effective theory is promoted by Gabaix et al Jit], based on the 
empirical power-law relation between the large trading volume and the large price return. Instead, from the phe- 
nomenological analysis of the order book data, Farmer and Weber et al. suggest that large price movements are driven 
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by the fluctuation in liquidity l32ll33ll . e.g., variations in the response to changes in supply and demand, and the low 
density of limit orders, etc. In fact, it is highly nontrivial to offer a complete answer, how large price movements occur. 
Nevertheless, it remains very important to fully understand the statistical properties of the financial fluctuations, such 
as the power-law distributions of different quantities, and their relations, as well as the long-range time correlations. 

In this paper we develop a multi-agent model of trading activity, aiming at a full understanding of the statistical 
properties of the financial fluctuations including the power-law distributions of the price return, trading volume and 
number of trades and the long-range time correlation of the volatility. In the literature, some stochastic models of 
trading activity, typically at the phenomenological level, have been analyzed for this purpose (34, 35, 3^]. Certain 
aspects of the financial fluctuations could be reproduced. A multiplicative stochastic model of the time interval 
between two successive trades, for example, is able to reproduce the statistical properties of the number of trades 
i 3411 . but does not refer to the relation with the power-law distributions of the price return and trading volume. In the 
present paper, we construct our model based on the microscopic structure and interactions, to capture fundamental 
mechanisms in the financial dynamics. 

The paper is organized as follows. In Sec. II, the dynamic herding model with interactions of trading volumes is 
introduced. In Sec. Ill, numerical results of the model are presented. Sec. IV contains the conclusion. 



2. Model definition 

The concept of percolating or herding is important in describing the financial markets 1 1^, 18, 37]. The dynamic 



version of the static percolation model, the so-called EZ herding model, shows certain attractive features 11811 . e.g., the 
herding structure is dynamically generated in a simple but robust way. The EZ herding model captures the power-law 
distribution of the price return, but the volatility is short-range correlated in time. To achieve the long-range time 



correlation of the volatility, a feed-back interaction should be introduced 1201I38H . Up to now, however, it is still far 
from realistic and harmonic. For example, only the volatility is concerned. The trading volume and number of trades 
have not been touched. The price return is calculated from the volatility with random ± 1 signs, and this should not 
describe the realistic financial dynamics. 

In this section, we develop a dynamic herding model including the price return, trading volume, and number of 
trades. In fact, it is not easy to build such a model. We have probed many possible variations of the microscopic 
structure and interaction, and finally come to the present form. 

2.1. Standard EZ herding model 

To start, let us first consider the EZ herding model. The system consists of M agents, which form clusters during 
dynamic evolution. Initially, each agent is a cluster. The dynamics evolves in the following way: 

(1) At a time step ?, select an agent i (and thus its cluster) at random. 

(2) With a probability 1 - a, i remains inactive in trading, and select another agent j randomly. If i and j are in 
different clusters, combine the two clusters into one. 

(3) With a probability a, i becomes active and makes a trade. Then all agents in the cluster follow. After that, this 
cluster is broken into a state that each agent is a separate cluster. The size of this cluster is recorded as s(t). 

Here the probability a is a constant, and controls the dynamic evolution. Since one does not define buying or 
selling of the trade, only the magnitude of the price return defined as |r(f)| = s(t) is essentially generated. The step 
(2) represents transmission of information. Considering the time between two actions as the time unit, I /a is the 
rate of transmission of information. If a is small, for example, transmission of information is fast, and agents tend to 



form larger clusters and act collectively. Numerical simulations [18] show that for a certain value of a, the probability 
distribution P(s) obeys a power law. 

2.2. Herding model interacted with trading volume 

However, the EZ herding model does not exhibit the long-range time correlation of the volatility. Furthermore, 
we need to include the price return, trading volume and number of trades (8, 23, 25]. Therefore, we assume that 



each agent trades with an individual trading volume v,(f) > 0. Denoting buying and selling with <x,(f) = +1 and -1 
respectively, all agents in a cluster are given a same trade sign <x,(f). Initially, each agent i with v,(0) = 1 is a cluster, 
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and randomly selects a trade sign cr,(0) — +1. The total trading volume is set to V(0) = 1. We construct the dynamics 
as following: 

(1) At a time step t, select an agent ; (and thus its cluster) at random, and calculate the trading probability 

alt) = , t' > 1. (1) 

l + bV(t- l)/vlt-f) 

Here v,(f - t') is the trading volume of ; at its last trade, V(t - 1) is the total trading volume at t - 1, and the parameter 
b is a positive value. 

(2) With a probability 1 - alt), the agent i remains inactive in trading, and select another agent j randomly. If i 
and j are in different clusters, combine the two clusters into one. The trade sign of the new cluster is taken to be that 
of the larger cluster of the previous two. 

(3) With a probability alt), all agents in the cluster which ; belongs to and another randomly selected cluster 
become active, and make trades according to their trade signs. After that, these two clusters are broken into a state 
that each agent is a separate cluster with a trade sign selected randomly. 

Our plausible observation is that if an agent is collecting much information, i.e., with a small alt), it may perform 
a large trade. Therefore, we assume v,(f) = 1 1 alt) to be the trading volume of the agent i. Then V(t) = 2/ vlt) is the 
total trading volume, with the sum over the two active clusters. Let N(t) denote the number of agents in the two active 
clusters, we define it as the number of trades. 

To determine the price return, we need more careful consideration. Empirical studies show that the trading volume 
seems to have a square root impact on the price return, and the price return saturates at extremely large trading volumes 
2311. Further, the correlation between the price return and trading volume is largely due to the number of trades 
39|] .Following the square root price impact function, we assume that the price return is determined by the volume 



imbalance Q(t) = Yi v,(f)o"/ and the number of trades. The volume imbalance reflects the difference between supply 



and demand 0321 13311 . Quantitatively, we define the price return 



lit) = S ign(Q(t)) l£l Ji - V#(0 ■ (2) 



The parameter A is taken to be a large positive value, such that r(f) ~ VICKOI y/N(t) at relatively small \Q(t)\, and 
r(t) ~ y/N(t) at extremely large \Q(t)\. 

The key ingredient in our model is the time-dependent probability alt). In Eq. ([1]), we assume that alt) depends 
on the ratio of the total trading volume at time t - 1 to its individual trading volume at its last trade. In financial 
markets, a large trading volume is usually accompanied by the strong fluctuation of the price return JI, 23, 25]. This 



inversely leads to large trading volumes in next time steps. Therefore, alt) is taken to be inversely proportional to 
the trading volume V{t - 1). If V{t - 1) is large, transmission of information is fast, the probability of combining 
two clusters is high, then the number of the trades increases on average, and finally leads to large trading volumes. 
Such a dynamic feed-back interaction of the trading volume essentially generates the long-range time correlation of 
the volatility. On the other hand, alt) is taken to be proportional to the individual trading volume vlt — f) at its last 
trade, based on the empirical assumption that an agent with a large trading volume in its last trade may be more active 
in trading in next time steps. Further, the thermodynamic limit is well defined due to the ratio V(t - l)/v,(f - f) in 
Eq. ©. 



3. Simulation results 

3.1. Probability distribution function 

In our model, the only tunable parameter is b. In calculating the price return, we fix the value A = 50. For each 
value of b, we take an average over 10 s iterations, after 10 6 iterations for equilibration. To detect the finite size effect, 
we perform extensive simulations with different total numbers of agents, and find that the results become stable for 
M > 40000 as shown in Fig. Q](a). Therefore, we report the results with M = 80000. 
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From empirical studies of the stock time series, the probability distribution of the trading volume V obeys a power 

law 

P(V) ~ y-( 1+ ^\ (3) 

with fa = 1.5, while that of the number of trades N obeys 

P(N) ~ AT( 1+ ?«\ (4) 
with fa — 3.4. The exponents of these two power-law distributions appear to have an approximate relation fa - 2 fa 

i. 

The probability distributions of the trading volume and number of trades in our model are carefully investigated. 
In Fig. [f|(a) and (b), P(V) and P(N) are plotted for b = 0.30, 0.45, 0.60 in log-log scale. For a small b, e.g., b < 0.3, 
P(V) and P (AO decay rapidly, and do not show a power-law behavior. As b increases, both P(V) and P(N) show a 
power-law behavior at b — 0.45, at least in two orders of magnitude. In this sense, the system exhibits a 'cross-over' 
behavior. Fitting the curves with the power laws in Eqs.(f3|l and ©, we estimate fa = 0.97 and fa = 2.1 1. These val- 
ues of the exponents are consistent with the approximate relation fa = 2fa. For b > 0.45, P(V) remains a power-law 
behavior in a broad range of b, but the exponent fa changes with b and becomes smaller. On the other hand, P(N) 
deviates from a power-law behavior up to a medium value N, while a power-law tail is still kept with an exponent fa 
approximately the same as that at b = 0.45. In Fig. [JJ the curves for b = 0.60 are displayed. 
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Figure 1: (Color online) Probability distributions of the trading volume and number of trades are plotted for b = 0.30, 0.45 and 0.60 in (a) and (b). 
The results are obtained with a total number of agents M = 80000. For comparison, a curve of b = 0.45 with M = 40000 is also displayed in (a). 

It is well known from the empirical analysis that the probability distribution of the price return exhibits a power- 
law tail 

PQr\) ~ k|- (1+ H (5) 

with g r = 3.0 10]. In Fig. [2] the probability distribution of the price return of our model is plotted for b = 0.45 in 
log-log scale. The curve can be nicely fitted with a power law, and the exponent is estimated to be 1.95, close to 
fa — 2.11. In summary, our dynamic herding model at b — 0.45 captures the power-law distributions of the trading 
volume, number of trades and price return. The corresponding exponents of these power-law distributions follow an 
approximate relation g r = fa = 2 fa. This is in agreement with that of the real markets reported in Ref. Q8Q, though 
the exponents themselves are slightly different. 

The specific form of a,-(f) is not very important. Other functions may also work, if their behaviors are similar to 
that in Eq. ([TJ. To verify this we also study the model with another a,(f) with a form 

fli(f) = 1 -ce~ d $%, t'>l, (6) 

where c and d are two positive parameters. We remain the inverse relation between a,(?) and the ratio V(t- l)/v,-(f — t'), 
thus makes it behave similar to Eq. ([TJ. An important characteristic of this model is the simple relation v,(f) = 1 /a,(f) 
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Figure 2: (Color online) Probability distribution of the volatility for b = 0.45 and M = 80000. 



between the trading volume and trading probability. It indicates that v,(0 are the dynamic variables, interact each 
other through the trading probability a,(f), and evolve according to Eq. ©■ 

We study the probability distributions of trading volume, number of trades and price return for the model with a,(f) 
defined in Eq. ([SJ. To make the probability a,(f) fixed between and 1, we set c = 1.0. By adjusting the parameter 
d, one observes a 'cross-over' behavior similar to the model with a,(f) defined in Eq. (|T): for small c, no power-law 
behavior is observed, while for large c power-law behavior occurs. In Fig. [3] power behaviors of P(V), P{N) and 
P(\r\) for a large c = 2.0 are plotted. By fitting the slopes of these curves, we estimate £y = 0.86, fjy = 1-89 and 
= 1.87. These exponents display a relation = g N = 2% v also consistent with that of the real markets. 




Figure 3: (Color online) Probability distributions of trading volume, number of trades, and price return for the model with Uj(t) defined in Eq. {6j. 
Results are obtained with c = 1.0, d = 2.0 and M = 80000. 



3.2. Auto-correlation function 

The long-range time correlation of the volatility is another important feature of the financial markets. Let us define 
the the auto-correlation function of the magnitude of the price return as 

<k(0||r(f + T)|>-(<|r(Q| » 2 

< KOI 2 > -(< KOI » 2 ' 

where < ■ ■ ■ > represents the average over the time t. From the empirical analysis, it obeys a power law, 

C(t) ~ T- 1 , (8) 



with ,1s 0.3 UG3]. 

We calculate C(t) of the model with a,(f) defined in Eq. ([]]), and observe that a power-law behavior is achieved 
only when b is in the very neighborhood of 0.45. In Fig. H](b), C(t) is plotted for b = 0.45 in log-log scale. The 
curve can be nicely fitted by a power law, indicating a long-range time correlation of the volatility. The exponent A is 
estimated to be 0.27, very close to that of the real markets. This improves the result A = 0.90 in a naive model 1 3811 . A 
power law with an exponent 0.53 (not shown in figure) is also observed in C(t) of the model with another a,(?) defined 
in Eq. ((6]), but with a reactively bigger fluctuation. 
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Figure 4: (Color online) auto-correlation function of the volatility for b = 0.45 and M = 80000. 

In our model, a large price movement is induced by a large volume imbalance \Q(t)\ and a large number of trades 
N(t) due to the price dynamics defined in Eq. (0. This may occur when large clusters exist in the system. Large 
clusters are formed after agents are active in collecting information for a period of time. However, a single large 
trading volume does not necessarily lead to a large price movement. 



4. Conclusion 

We introduce a dynamic herding model with interactions of trading volumes. At time t, each agent trades with a 
probability ai(t) as a linear function of the ratio of the total trading volume V(t - 1) at time t - 1 to its own trading 
volume Vi(t - f') at its last trade. Agents are endowed with trade signs cr,(f) = +1 denoting buying and selling, and 
the price return is determined by the volume imbalance and number of trades. We find that at a threshold b = 0.45, 
our model reproduces the power-law distributions of the trading volume, number of trades and price return, and more 
importantly, the approximate relation f r = fa = 2£y. The generated volatilities are long-range correlated in time. 
We also investigate the model with another exponential form of a,(f) which remains its inverse relation with the 
ratio V(t - l)/v,-(f - f')> and similar power-law behaviors and long-range correlation are observed. This indicates the 
robustness of our model. 

By simulating the agents' reactions to the market information through the interactions between their individual 
trading volumes and the total trading volumes, the model can explain the power-law behaviors of the trading volume, 
number of trades and price return, and their relations. We believe that our model captures certain essences of the 
financial markets. Further works should include understanding of the Leverage and anti-Leverage effects in western 
and Chinese financial markets |4(| 41, 42, 43], and more general interactions between the price return, trading volume 



and number of trades, etc. 
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